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Evaluation of endohedral doping of hydrogenated Si fullerenes as a route to magnetic 

Si building blocks 
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Density-functional theory based global geometry optimization is used to scrutinize the possibility 
of endohedral doping of hydrogenated Si fullerenes as a route to Si nanostructures with high magnetic 
moments. In contrast to previous suggestions, our unbiased sampling finds the smallest Sii6Hi6 
endohedral cage generally too small to encapsulate 3d metal dopant atoms. For the next larger 
fullerene-like cage though, we identify perfectly symmetric A/@Si2oH2o (M=Co, Ti, V, Cr) cage 
structures as ground states. These structures conserve the high spin moment of the dopant atom 
and therewith underscore the potential of this Si nanoform for novel cluster-based materials with 
unique magnetic properties. 

PACS numbers: 61.46.Bc, 36.40.Cg, 73.22.-f, 31.15.es 



The prevalence of silicon in semiconductor industry has 
raised great interest in silicon-based nanostructures that 
could act as building blocks for new easy-to- integrate and 
engineered materials. In contrast to the compact and 
reactive form of pure Si clusters, endohedral metal dop- 
ing was found to produce saturated fullcrcnc and other 
polyhedral cage structures, and therewith an intriguing 
novel nanoform of Si with great promise e.g. for opto- 
electronic device applications [1-14]. In these structures 
the cage geometry is stabilized through a strong interac- 
tion with the encapsulated dopant atom, which unfortu- 
nately goes hand in hand with a quenching of the dopant 
spin moment [15]. The recent suggestion that also hy- 
drogen termination of Si clusters could yield an (empty) 
fullcrcnc configuration [16-18] has thus raised hopes 
that metal-doping corresponding hydrogenated fullerenes 
would yield cage structures with minimized M-Si inter- 
action [18]. With the atomic character of e.g. magnetic 
dopants then likely conserved, this would offer a route to 
specifically develop Si fullerene species with large mag- 
netic moments for magneto-electronic applications. 

In order to further explore and substantiate this idea 
we perform density-functional theory (DFT) based global 
geometry optimizations of the smallest hydrogenated en- 
dohedral cage, Sii6Hi6, for a range of metal dopants. In 
order to address a possible influence of the size of the 
dopant atoms, this range comprises not only previously 
discussed high magnetic moment species like Ti or Cr 
[18], but extends over the entire 3d series. While our ex- 
tensive unbiased sampling indeed confirms the cage-like 
geometry as global ground state of empty Sii6Hi6, this 
is unfortunately not the case for any of the 3d dopants. 
Instead of a symmetric metal encapsulation, strongly dis- 
torted or even broken cages are significantly more stable 
in all cases. Moreover, the stronger M-Si interaction in 
these ground-state structures leads to the same quench- 
ing of the magnetic moment as in the regular M@Sii6 
fullerenes. Attributing these discouraging findings to the 
insufficient space inside Sii6Hi6, we therefore extend our 



study to the next larger fullerene-like structure, Si2oH2o 
[16]. Testing the high magnetic moment sequence Ti, 
V, Cr, and Co we now obtain in all cases perfectly sym- 
metric endohedral cage geometries as ground states. Spin 
density distribution analysis furthermore reveals that the 
high spin states of these structures come indeed primarily 
from conservation of the dopant magnetic moment, i.e. 
the tested sequence already offers an intriguing toolbox 
of Si nanoforms with septet to quartet spin, respectively. 
Confirming also the next larger Si24H24, Si26H26 and 
Si28H28 cages in the fullerene sequence as correspond- 
ing ground-state geometries, this strongly supports the 
idea of generating promising magnetic building blocks 
by doping sufficiently-sized hydrogenated Si cages. 

All ground-state total energy calculations in this work 
have been performed with the all-electron full-potential 
DFT code FHI-aims [19], treating electronic exchange 
and correlation within the generalized-gradient approxi- 
mation functional due to Pcrdcw, Burke and Ernzerhof 
(PBE) [20] . All sampling calculations were done with the 
" tier2" atom-centered basis set and using " tight" settings 
for the numerical integrations. Stability of the identified 
minima has been confirmed by frequency analysis. For 
the ensuing electronic structure analysis of the optimized 
geometries, the electron density was recomputed with an 
enlarged "tier3" basis set [19]. As in our previous work on 
metal-doped Si clusters [13, 15] systematic convergence 
tests indicate that these settings are fully converged with 
respect to the target quantities, here the energetic differ- 
ence of isomers, as well as total electron and spin density 
distributions. For the actual global geometry optimiza- 
tion procedure we relied on the basin-hopping (BH) al- 
gorithm [21, 22], which samples the huge configurational 
space through consecutive jumps from one local mini- 
mum of the potential energy surface to another. In our 
specific BH implementation [13, 15, 23] we achieve this 
by randomly moving one or all atoms in the cluster to 
a new position within a sphere of 1.5 A, each time fol- 
lowed by a local geometry optimization and acceptance 
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or rejection of the thus created structure on the basis of 
a Metropolis algorithm. 

The starting point of our investigation is an extended 
configurational search, which even after thousands of 
trial structures confirms that the ground-state structure 
of Sii6Hi6 is indeed the perfectly symmetrical tetrahcdral 
(Trf) cage shown in Fig. 1 [16]. Featuring Si atoms at 
distances of 2.95 A and 3.17 A from the center, this cage 
provides in principle enough geometric space to host an 
cndohcdral dopant atom, with the hydrogenation hoped 
to minimize the actual M-Si interaction. Unfortunately, 
the global geometry optimization results for MSii6Hi6 
clusters with metal dopant atoms ranging over the en- 
tire 3d series from Sc to Zn reveal a different picture. 
In all cases only a few BH steps are required to identify 
significantly more stable geometries, which correspond 
to severely distorted cages or more compact structures 
that push the dopant atom from the center to the cluster 
fringe (or even expel it completely in case of Zn), cf. Fig. 
1 [24]. With the exception of Zn (0.37 eV) the energetic 
difference of these new structures to the symmetric endo- 
hedral cage is of the order of 1 eV or larger, which makes 
a kinetic trapping in the high-energy tetrahedral cage ge- 
ometry rather unlikely. It also suggests that the obtained 
energetic order is not affected by the well-known limi- 
tations of the employed semi-local DFT functional. To 
further validate this, we nevertheless evaluated the en- 
ergetic differences also with the hybrid PBEO functional 
[25], without obtaining any qualitative changes. 

As already indicated by the typically high Si coordi- 
nation, the identified low-energy isomers exhibit a rather 
strong M-Si interaction involving substantial charge rear- 
rangements. Not surprisingly, this goes along with an al- 
most complete quench of the dopant spin moment. Clus- 
ters having an even number of electrons are singlets (Ti, 
Cr, Fe, Ni and Zn doping), clusters with odd number of 
electrons are doublets (Sc, V, Mn, Co and Cu doping), 
with again no changes observed when repeating the calcu- 
lations at the PBEO hybrid functional level. In case of the 
doublets, the unpaired electron is hereby not even neces- 
sarily centered on the dopant atom: For both Co and Sc 
the spin density is either partly or even mainly located on 
Si atoms, respectively Obviously, this completely rules 
out the original proposition of using these nanomatcrials 
as high magnetic moment Si building blocks. 

Suspecting an insufficient space inside the cage as rea- 
son for these findings we proceed to the next larger hy- 
drogenatcd fullerene-like cage Si2oH 2 o, for which exten- 
sive sampling again confirms the empty cage geometry 
shown in Fig. 2 as ground-state structure. Intriguingly, 
the increased radius of 3.33 A of this icosahedral (Ih) cage 
now seems to be large enough to even accommodate most 
space-demanding high spin-state dopant atoms. Focus- 
ing on the sequence Co, Ti, V, and Cr with quartet, quin- 
tet, sextet and septet atomic spin moments, respectively, 
the cndohcdrally-dopcd symmetric cage corresponds each 



time to the most stable structure found in our extensive 
sampling runs. Moreover, the energetic gap to the second 
lowest identified isomer, which then corresponds to a dis- 
torted or broken cage structure, is in all cases larger than 

0. 5 eV. From this we again do not suspect our findings to 
be jeopardized by the limitations of the employed semi- 
local DFT functional; an interpretation that we find con- 
firmed by hybrid PBEO-level calculations that yield for 
all four dopant atoms the same energetic order of ground- 
state cage and second lowest isomer as in the PBE case. 

For this fullercne size we thus find a beautiful confirma- 
tion of the hydrogenation concept, considering that previ- 
ous global geometry optimization work clearly identified 
the corresponding endohedral M@Si 2 o cages (M =Ti, V, 
Cr) as metastable, with the ground-states of these un- 
hydrogenated clusters instead given by heavily distorted 
and reactive structures [13]. In contrast to the again al- 
most complete spin quench in the latter structures, the 
doped hydrogenated fullerencs furthermore perfectly con- 
serve the high magnetic moments of their metal dopants, 

1. e. Co@Si 20 H 20 , Ti@Si 20 H 20 , V@Si 20 H 20 , Cr@Si 20 H 20 
exhibit quartet, quintet, sextet and septet spin mo- 
ments, respectively. Analysis of the spin density distri- 
bution confirms that the unpaired electrons indeed re- 
side predominantly on the central metal atom, which is 
in line with the obtained almost negligible hybridization 
of the occupied metal 3d states. The metal-cage inter- 
action instead seems largely mediated via the metal 4s 
state, which rationalizes the stabilization of a symmetric 
dopant position in the center of the cage. 

Despite considerable binding energies, which at the 
PBE-level exceed 1 eV [24] and therewith stabilize the 
cndohcdral cage geometry, the atomic character of the 
metal dopant atom is thus largely conserved. This is 
also apparent in the optical excitation spectra, computed 
within the time-dependent DFT (TD-DFT) linear re- 
sponse formalism as implemented in the Gaussian03 suite 
[26, 27]. Figure 3 illustrates this for the Cr@Si 20 H 20 case. 
In stark contrast to corresponding spectra of e.g. the bro- 
ken cage structures of MSii6Hi6, only a few well-defined 
transitions arc apparent in the optical range. These tran- 
sitions are energetically shifted and changed in their or- 
der in comparison to the corresponding spectra of the 
isolated atomic dopants. Nevertheless the "atomic" char- 
acter is preserved, i.e. apart from their magnetic proper- 
ties the endohedrally-doped hydrogenated fullerenes also 
exhibit promising optical properties. 

In summary, we have thus systematically scruti- 
nized the possibility of stabilizing endohedrally doped Si 
fullerenes through hydrogenation. Our unbiased global 
geometry optimization based on semi-local and hybrid 
functional DFT undoubtedly revealed that the smallest 
hydrogenated Sii6Hi6 cage focused on in previous work 
is generally too small to encapsulate 3d transition metals. 
On a conceptual level this underscores the danger of bas- 
ing nanoscalc materials design through predictive-quality 
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FIG. 1: Identified ground-state geometry of Sii6Hi6, as well as of M@Sii6Hi6 clusters (with M covering the 3d metal series) 
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FIG. 2: Identified hollow cage ground-state geometries of the 
next larger hydrogenated Si fullerenes. 



theory on incomplete explorations of the vast configura- 
tional space, like simple comparisons of chemically intu- 
itive candidate structures. For the next larger fullerene- 
like cage though, our first-principles sampling indeed 
identifies perfectly symmetric M@Si2oH 2 n (M=Co, Ti, 
V, Cr) cage structures as ground states. Compared to 
distorted or broken cage geometries, these structures are 
stabilized through substantial delocalized 4s-cage inter- 
action, while nevertheless largely retaining the atomic 
character of the metal dopant. With respect to mag- 
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FIG. 3: Calculated TD-DFT optical excitation spectrum of 
the endohedral Cr@Si2oH2o fullerene ground-state, with the 
insets illustrating the dominant Kohn-Sham orbitals behind 
the individual transitions. 



netic properties, the confirmed quartet, quintet, sextet 
and septet spin moments of Co@Si2oH 2 o, Ti@Si2oH 2 o, 
V@Si2oH2o, and Cr@Si 2 nH 2 n, respectively, thus already 
offer a nice toolbox of unreactive building blocks with 
high magnetic moments. Nicely shielded by the fullerene 
cage, the encapsulated dopants furthermore also conserve 
atom-like optical properties. This strongly suggests en- 
dohedral doping of hydrogenated fullerenes as a viable 
route to novel cluster-based materials for magneto-optic 
applications. This vision gets further support by exten- 
sive sampling runs that also suggest empty cage geome- 
tries as ground-states for the next larger hydrogenated 
Si fullerenes up to Si28H 2 g, cf. Fig. 2. Eventually, 
these cages will offer enough space to also host multi-core 
dopants. Encapsulation into hydrogenated Si fullerenes 
thus appears as a promising avenue to generally protect 
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metal clusters from reactive environments and make their 
unique material's properties available in applications. 
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